(MT@ DDTI original research article 



MICROBIOLOGY 



published: 27 December 2013 
doi: 10.3389/fmicb.2013.00414 




Updating the Vibrio clades defined by multilocus sequence 
phylogeny: proposal of eight new clades, and the 
description of Vibrio tritonius sp. nov. 

Tomoo Sawabe 1 *, Yoshitoshi Ogura 2 , Yuta Matsumura 1 , Gao Feng 1 , AKM Rohul Amin 1 , 
Sayaka Mino \ Satoshi Nakagawa 1 , Toko Sawabe 3 , Ramesh Kumar 4 , Yohei Fukui 5 , 
Masataka Satomi 5 , Ryoji Matsushima 5 , Fabiano L. Thompson 6 , Bruno Gomez-Gil 7 , 
Richard Christen 8 - 9 , Fumito Maruyama 10 , Ken Kurokawa 11 and Tetsuya Hayashi 2 

' Laboratory of Microbiology, Faculty of Fisheries Sciences, Hokkaido University, Hakodate, Japan 

2 Division of Genomics and Bioenvironmental Science, Frontier Science Research Center, University of Miyazaki, Miyazaki, Japan 

3 Department of Food and Nutrition, Hakodate Junior College, Hakodate, Japan 

4 National Institute for Interdisciplinary Science and Technology (CSIR), Kerala, India 

5 National Research Institute of Fisheries Science, Fisheries Research Agency, Yokohama, Japan 

6 Department of Genetics, Center of Health Sciences, Federal University of Rio de Janeiro (UFRS), Rio de Janeiro, Brazil 

7 A.C. Unidad Mazatlan, CI AD, Mazatlan, Mexico 

8 CNRS UMR 7138, Systematique-Adaptation-Evolution, Nice, France 

9 Systematique-Adaptation-Evolution, Universite de Nice-Sophia Antipolis, Nice, France 

'" Graduate School of Medical and Dental Sciences, Tokyo Medical and Dental University, Tokyo, Japan 
" Earth-Life Science Institute, Tokyo Institute of Technology, Tokyo, Japan 



Edited by: 

Daniela Ceccarelli, University of 
Maryland, USA 

Reviewed by: 

Xiu-Lan Chen, Shandong University, 
China 

Matteo Spagnoletti, University 
College London, UK 

'Correspondence: 

Tomoo Sawabe, Laboratory of 
Microbiology, Faculty of Fisheries 
Sciences, Hokkaido University, 
3-1-1 Minato-cho, Hakodate 
041-8611, Japan 

e-mail: sawabe@fish.hokudai.ac.jp 



To date 142 species have been described in the Vibrionaceae family of bacteria, classified 
into seven genera; Aliivibrio, Echinimonas, Enterovibrio, Grimontia, Photobacterium, 
Salinivibrio and Vibrio. As vibrios are widespread in marine environments and show 
versatile metabolisms and ecologies, these bacteria are recognized as one of the most 
diverse and important marine heterotrophic bacterial groups for elucidating the correlation 
between genome evolution and ecological adaptation. However, on the basis of 16S 
rRNA gene phylogeny, we could not find any robust monophyletic lineages in any of 
the known genera. We needed further attempts to reconstruct their evolutionary history 
based on multilocus sequence analysis (MLSA) and/or genome wide taxonomy of all the 
recognized species groups. In our previous report in 2007, we conducted the first broad 
multilocus sequence analysis (MLSA) to infer the evolutionary history of vibrios using nine 
housekeeping genes (the 16S rRNA gene, gapA, gyrB, ftsZ, mreB, pyrH, recA, rpoA, and 
topA), and we proposed 14 distinct clades in 58 species of Vibrionaceae. Due to the 
difficulty of designing universal primers that can amplify the genes for MLSA in every 
Vibrionaceae species, some clades had yet to be defined. In this study, we present a better 
picture of an updated molecular phylogeny for 86 described vibrio species and 10 genome 
sequenced Vibrionaceae strains, using 8 housekeeping gene sequences. This new study 
places special emphasis on (1) eight newly identified clades (Damselae, Mediterranei, 
Pectenicida, Phosphoreum, Profundum, Porteresiae, Rosenbergii, and Rumoiensis); (2) 
clades amended since the 2007 proposal with recently described new species; (3) 
orphan clades of genomospecies F6 and F10; (4) phylogenetic positions defined in 3 
genome-sequenced strains (N418, EX25, and EJY3); and (5) description of V. tritonius sp. 
nov., which is a member of the "Porteresiae" clade. 



Keywords: vibrios, Vibrionaceae, multilocus sequence analysis, evolution, housekeeping protein gene. Vibrio 
tritonius 



INTRODUCTION 

Bacterial systematics has evolved alongside the development of 
innovative methodologies and techniques (Wayne et al., 1987; 



Abbreviations: AFLP, amplified fragment length polymorphism; ANI, average 
nucleotide sequence identity; AAI, average amino acid identity; DDH, DNA-DNA 
hybridization; MLSA, multilocus sequence analysis; MLST, multilocus sequence 
typing; MP, maximum parsimony; NJ, neighbor joining; ML, maximum likelihood; 
PFGE, pulse field gel electrophoresis; RAPD, random amplified polymorphic DNA. 



Stackebrandt et al., 2002; Gevers et al., 2005). The first definition 
of bacterial species in "phylogenetic terms" was developed in 1987 
using the DNA-DNA reassociation and DNA sequencing. These 
approaches to bacterial systematics provided us with a uniform 
definition of prokaryotic species (Wayne et al., 1987). In 2002, an 
ad hoc committee listed additional innovative methods that could 
be used for bacterial systematics, such as 16S rRNA gene sequence 
analysis, DNA typing methods (AFLP, RAPD, Rep-PCR, PFGE), 
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MLSA, WGS analysis, and proteomics (Stackebrandt et al., 2002). 
The primary purpose of the committee's statement was to pro- 
mote dialogue among systematists, population and evolutionary 
geneticists, ecologists and microbiologists for the benefit of bac- 
terial systematics in general, and to create a more transparent 
species concept in particular (Stackebrandt et al., 2002). Among 
those innovative methodologies, MLSA and the WGS analysis 
have become the most important and successful methodologies; 
their strong impact on bacterial systematics is due to data repro- 
ducibility and portability (see, e.g., Maiden et al., 1998; Aanensen 
and Spratt, 2005; Gevers et al., 2005; Konstantinidis and Tiedje, 
2005; Staley, 2006; Goris et al, 2007; Richer and Rossello-Mora, 
2009; Auch et al, 2010). 

MLST, the prototype for MLSA-based methodology, was used 
for the first highly portable typing of Neisseria meningitides from 
invasive disease and healthy carriers, and it yielded the first under- 
standing of the epidemiology and population structure of that 
infectious agent (Maiden et al., 1998). Its high levels of discrimi- 
natory power between those strains, which required half the loci 
typically required for a classical allozyme electrophoresis, and its 
superior application to evolutionary, phylogenetic, or population 
genetic studies, allowed researchers to develop MLST schemes 
for a number of bacteria taxa (Aanensen and Spratt, 2005) (also 
refer to the MLST website; http://www.mlst.net/). It also opened 
the use of MLSA for bacterial systematics (e.g., Sawabe et al., 
2007; Thompson et al, 2007; Bishop et al, 2009), and it guided 
the reconsideration and re-evaluation of prokaryotic species con- 
cepts (Gevers et al., 2005; Staley, 2006; Preheim et al, 2011). 
Even now, MLSA provides a better understanding of taxonomi- 
cally controversial bacterial taxa; for example, the human origins 
of the Agrobacterium (Rhizobium) radiobacter clustered as a well- 
separated genovar (Aujoulat et al., 2011) and the highly versatile 
aeromonads consisting of 3 major clades (Roger et al., 2012). 

In the genome era, genome sequencing has been used to 
characterize new bacterial species (Haley et al, 2010; Hoffmann 
et al., 2012), to reclassify bacterial taxons such as Neisseria 
(Bennett et al, 2012), Acinetobacter (Chan et al., 2012) and Vibrio 
(Lin et al., 2010), and to challenge defined prokaryotic species 
(Konstantinidis and Tiedje, 2005; Thompson et al., 2009; Chan 
et al., 2012). Using these WGSs, in-silico DDH calculations can 
also be emulated, mainly in two ways: high-scoring segment pairs 
(HSPs) (Konstantinidis and Tiedje, 2005; Goris et al, 2007) and 
the genome-to-genome distance calculation, called the Digital 
DDH measurement (Auch et al., 2010). The criterion of more 
than 95% ANI is currently a widely used similarity value for 
species delineation. The WGS analysis also supercedes the limita- 
tions of MLSA, which is only capable of including genes that are 
successfully amplified by designed primers (Gevers et al, 2005; 
Thompson et al, 2005; Sawabe et al, 2007). 

Vibrionaceae are at the forefront of bacterial taxons being 
tested with new innovative methodologies and techniques for 
bacterial systematics (Thompson et al, 2001, 2005, 2009; Sawabe 
et al., 2007). The number of species described in Vibrionaceae 
has increased remarkably since the establishment of genome 
fingerprinting techniques (Thompson et al., 2001) and MLSA 
schemes (Thompson et al, 2005, 2007; Sawabe et al., 2007, 
2009). Now, a total of 142 species are recognized in the family 



Vibrionaceae (Association of Vibrio Biologists website; http:// 
www.vibriobiology.net/). Vibrionaceae are defined as a group of 
strains with the following characteristics: they are Gram-negative 
rods with a polar flagellum enclosed in a sheath, have faculta- 
tive anaerobic metabolisms, are capable of fermenting D-glucose, 
and grow at 20° C. The bacteria are primarily aquatic, and most 
species are oxidase positive, can reduce nitrate to nitrite, require 
Na + for growth, and ferment D-fructose, maltose, and glycerol 
(Gomez-Gil et al., 2014). In addition, most vibrio species fer- 
ment a variety of carbohydrates without gas production, and 
grow on TCBS medium (Farmer et al., 2005; Thompson et al., 
2009). As we experience a rapid expansion in the number of 
known species in the family Vibrionaceae, we face a number of 
unique vibrio isolates that lack one or more of the above com- 
mon properties. Vibrionaceae species are metabolically versatile, 
and the number of species showing gas production, nitrogen fix- 
ation, phototrophy, and non-motility is increasing (Gomez-Gil 
et al., 2014). Considerably more attention should be paid to the 
biological and genetic plasticity of vibrios to help understand the 
dynamics of vibrio evolution (Sawabe et al, 2007; Grimes et al., 
2009; Thompson et al., 2009). 

Due to the limitations of using 16S rRNA gene phylogeny 
(Gomez-Gil et al., 2014) to elucidate an "Integrated Vibrio 
Biology" that includes a biodiversity assessment, an inferred 
evolutionary history, the population biology, and genomics, 
Sawabe et al. (2007) developed an MLSA scheme for the family 
Vibrionaceae using nine gene sequences (ftsZ, gapA, gyrB, mreB, 
pyrH, recA, rpoA, topA and the 16S rRNA gene). The analysis 
involved the complete sequence sets of 9 genes from 58 vibrio 
taxa, and it revealed 14 monophyletic clades with a significant 
bootstrap support. The species within each clade shared >20% 
DDH, <5% G+C (mol%), >85% MLSA sequence similarity, 
and >89% AAI (Sawabe et al, 2007). 

Recent extraordinary progress in biodiversity studies and WGS 
projects in vibrios has resulted a substantial leap in novel vibrio 
species and taxonomically unassigned strains. In fact, after the 
proposal of 14 robust clades in 2007, more than 60 new species 
have been described. The number of genome-sequenced strains 
has exceeded 1000. It is, therefore, obvious that the phylogenetic 
tree based on the multilocus gene sequences reported in 2007 
is insufficient to show the most recent molecular phylogenetic 
structure of vibrios. 

In our previous analysis of multilocus gene phylogeny, eight 
housekeeping protein coding genes and the 16S rRNA gene were 
included in the MLSA. However, the 16S rRNA gene has a rather 
low interspecies resolution (100, 99.5, and 95.3% of maximum, 
median, and minimum resolution, respectively) (see Figure 5S, in 
Sawabe et al., 2007). It was also difficult to include the 16S rRNA 
gene sequences for the calculation of the radiation time of each 
clade. We have been debating whether the inclusion of the 16S 
rRNA gene sequences is necessary to infer the evolutional history 
of vibrios, but there are no other proper and fast tools to check 
the tree topologies constructed using a rather large data set (58 
species). 

A variety of methods have been proposed to tackle the prob- 
lem of gene tree reconciliation to reconstruct a species tree. When 
the taxa in all the trees are identical, the problem can be stated 
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as a consensus tree problem (Guenoche, 2013). The comparison 
of gene trees and their assembly into a unique tree representing 
the species tree is a general problem in phylogeny. However, in 
this study, we faced a different problem: determining whether 
the inclusion or exclusion of a given gene in the analysis would 
substantially change the outcome. This was the method used in 
this analysis to investigate if the inclusion or exclusion of the 16S 
rRNA gene sequences in the M LSA analysis would, or would not, 
affect the final result. 

The aims of this study were to re-evaluate how the 16S rRNA 
gene sequence affects the final phylogenetic tree, as based on 
multilocus gene sequencing analysis; to update our knowledge in 
vibrio biodiversity and evolution on the basis of an 8-gene MLSA; 
and to reconstruct a better vibrio phylogeny. The analysis pro- 
vided a further opportunity to propose additional eight clades to 
the most up-to-date vibrio phylogeny. 

MATERIALS AND METHODS 

SUBTREE INCONGRUENCE TEST OF MULTILOCUS GENE PHYLOGENY 

The usual approach to compare trees is to count how many sub- 
trees they share; a given subtree often has a different topology 
according to the method used (NJ, ML, or MP) but the differences 
are often subtle and generally not well-supported. Accordingly, it 
is best to compare trees according to bipartitions (a tree is consid- 
ered as a set of bipartitions, each one corresponding to an internal 
edge of the tree, the external ones connecting the leaves to the 
tree) (Guenoche, 2013). This was the method used in this anal- 
ysis to investigate whether the inclusion or exclusion of the 16S 
rRNA gene sequences in the MLSA analysis would, or would not, 
affect the final result. TreeDyn (Chevenet et al, 2006) was used to 
compare trees to subtrees sharing the same topology. A dedicated 
python script (using libraries from Huerta-Cepas et al, 2010) was 
used to compare trees with shared sub-trees, independently of 
their topologies. 

Our previous MLSA revealed that the 16S rRNA gene, in 
contrast to the other gene sequences used, has a rather low inter- 
species resolution (Sawabe et al., 2007). To increase the sensitivity 
and reduce the time taken by MLSA to update the vibrio phy- 
logeny, we compared the subtree topologies obtained from a 
nine-gene data set (that included the 16S rRNA gene) and an 
eight-gene data set (that excluded the 16S rRNA gene). We used 
the gene sequence data set from the 58 species, for which we had 
the sequence of every gene, and which is identical to the data set 
used in Sawabe et al. (2007). The method selected, (1) shared sub- 
trees only if they had the same topology, and (2) subtrees that 
shared the same species, independently of their topology. The 
results were visualized using TreeDyn (Chevenet et al, 2006), and 
congruent subtrees with the same topologies were indicated using 
the same color. 

SEQUENCING OF HOUSEKEEPING PROTEIN-CODING GENES 

An additional eight housekeeping genes of type strains of the 
genera Vibrio and Photobacterium were sequenced manually 
according to Sawabe et al. (2007) with newly designed primers 
( Vgap 1 50f: ACTCAYGGYCGTTTCAACGGYAC, Vgap957r:RC 
CGATTTCGTTRTCGTACCAAG, VftsZf55f:GTKGGTGGCG 
GCGGCGGTAA, VftsZ782r:ACACCACGWGCACCAGCAA 



GATCG, VftsZ-9f:ACCGATGATGGAAATGTCTGACGATG 
C, VmreB225f:RATGAAA GACGGCGTWATYGC, VmreB102 
5r:TCGCCRCCGTGCATRTCGATCA) (Table SI). All strains 
were maintained in ZoBell 2216E agar and stored with 20% 
glycerol at — 80°C. Whole genome sequencing was performed 
in nine type strains (V^ aerogenes LMG 19650 T , V. gazogenes 
ATCC 29988 T , V. halioticoli IAM 14596 T , V. neonatus HDD3-1 T , 
V. porteresiae MSSRF 30 T , V. rhizosphaerae MSSRF 3 T , V. ruber 
LMG 23124 T , V. tritonius sp. nov. AM2 T , and V. superstes 
G3-29 T ) using the Roche 454 FLX titanium genome sequencer 
alone or in combination with the Illumina MiSeq sequencer. 
The sequence reads were assembled using the Newbler soft- 
ware version 2.3 or later. Illumina reads were used only for 
sequence error correction. After auto-annotation by Microbial 
Genome Annotation Pipeline (MiGAP, http://www.migap.org/; 
Sugawara et al., 2009), relevant housekeeping gene sequences 
were retrieved and used for the MLSA. The house keeping genes 
necessary for updating the vibrio phylogeny were also retrieved 
from the latest version of the NCBI microbial genome and 
GenBank database (Release 197.0, 15 August 2013), and used 
in the analysis. All sequence data used in this study are listed in 
Table SI. 

SEQUENCE ANALYSIS 

MLSA was performed the same way as in Sawabe et al. (2007). 
The sequences were aligned using the ClustalX program (Larkin 
et al., 2007). The domains used to construct the phylogenetic trees 
shown in Figures 2, 3 were regions of the ftsZ, gap A, gyrB, mreB, 
pyrH, recA, rpoA, and topA genes of Vibrionaceae: positions 195— 
630, 225-861, 441-1026, 390-897, 171-543, 429-915, 87-873, 
and 570-990 (V. cholerae Ol Eltor N16961 (AE003852) num- 
bering), respectively. The regions were within those used in the 
previous study (Sawabe et al, 2007). Sequence similarity and the 
number of nucleotide and amino acid mutation were determined 
using MEGA version 5 (Tamura et al., 201 1). 

Split Decomposition Analysis (SDA) was also performed 
as described in Sawabe et al. (2007) using SplitsTree version 
4.12.8, with a neighbor net drawing and Jukes-Cantor correction 
(Bandelt and Dress, 1992; Huson and Bryant, 2005). The con- 
catenated sequences of the eight housekeeping genes were also 
generated using the program and used for a phylogenetic analysis 
combined with NJ, MP, and ML analyses (Sawabe et al, 1998). 

PHYLOGENETIC, GENETIC, PHEN0TYPIC, AND CHEM0TAX0N0MIC 
CHARACTERIZATION OF Vibrio tritonius sp. nov 

Four isolates of V. tritonius sp. nov., JCM 16456 T = LMG 
25401 T = AM2 T , JCM 16457 = LMG 25402 = MA12, JCM 
16458 = LMG 25403 = MAI 7, and JCM 16459 = LMG 25404 = 
MA35, isolated from gut of a sea hare, Aplysia kurodai, were used 
in this study. The strains were cultured on ZoBell 2216E agar 
(Oppenheimer and ZoBell, 1952) and stored at -80°C in 10% 
glycerol-supplemented broth. 

A total of 1400 bp 16S rRNA gene sequences of the four 
strains were determined according to Sawabe et al. (1998) using 
four sequence primers (24F, HOOF, 920R, and 1509R). The 16S 
rRNA gene sequences were blasted to the latest release ver. 197 
of GenBank and related sequences were retrieved. Finally, 16S 
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rRNA gene sequences of V. aerogenes X74705, V. brasiliensis 
AJ316172, V. cholerae X76337, V. fluvialis X74703, V. furnissii 
X76336, V. gazogenes X74705, V. hepatarius AJ345063, V. nereis 
X74716, V. porteresiae EF488079, V. rhizosphaerae DQ847123, 
V. ruber AF462458, V. tubiashiiX74725, and V. xuii AJ316181 were 
included in the phylogenetic analysis (Figure 4). Phylogenetic 
trees were constructed using three different methods (NJ, ML, 
and MP). For NJ analysis, distance matrices were calculated using 
the Kimura two parameters correction and using MEGA version 
5.0 (Tamura et al., 2011). ML and MP analysis was conducted 
using PHYLIP (Phylogeny Inference Package, version 3.573c, dis- 
tributed by J. Felsenstein, Department of Genetics, UW, Seattle, 
WA, USA). Sequences corresponding to positions 86-1420 of the 
E. coli gene (NC_000913) were used in this analysis. Figure 4 rep- 
resents a subset of the final tree obtained using the NJ method 
with 500 bootstrap replications. Nodes supported by ML and MP 
analyses are indicated by the bootstrap values in Figure 4. 

DNAs of bacterial strains were prepared following the proce- 
dures of Marmur (1961), with minor modifications. The mol% 
G+C content of DNAs was determined by HPLC (Tamaoka and 
Komagata, 1984). DNA-DNA hybridization experiments were 
performed in microdilution wells using a fluorometric direct 
binding method described by Ezaki et al. (1988); Ezaki et al. 
(1989). DNA-DNA similarity data were shown as the average 



value of triplicate experiments. V. brasiliensis, V. furnissii, V. flu- 
vialis, V. tubiashii, V. hepatarius, V. mytili, V. nereis, V. porteresiae, 
and V. xuii were selected as the reference species of these DNA- 
DNA hybridization experiments, based on the results from the 
MLSA molecular phylogenetic assessment and the 16S rRNA gene 
phylogeny of V. tritonius sp. nov. 

A total of 62 phenotypic characteristics were determined using 
the standard manual characterization method established in our 
laboratory (Sawabe et al., 1998). The carbon assimilation test 
was conducted using a basal seawater medium, as previously 
described (Sawabe et al, 1998). These phenotypic characteriza- 
tions were performed at 25°C. 0/129 sensitivity was determined 
using the sensitivity basal agar medium (Nissui Pharmaceutical, 
Tokyo, Japan) at 30°C. 

RESULTS 

THE SUBTREE INCONGRUENCE TEST BETWEEN THE 9-GENE AND 
8-GENE PHYL0GENIES 

Subtrees obtained in the 8-gene phylogeny were compared to 
those in the 9-gene phylogeny and the 16S rRNA gene phylogeny 
(Figure 1). Among the 13 subtrees reconstructed in the 8-gene 
phylogeny, 12 were retained in the 9-gene phylogeny; the only 
difference observed was the inclusion of V. proteolyticus in the 
V. cholera subtree (Figures 1A,B). Most of the 13 subtrees in the 
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V. proteolyticus LMG 3772 T 
V. harveyi LMG 4044' 
V. campbellii IMG 1 1 21 6 T 
V. rotiferianus LMG 21460 T 
V.natriegens LMG 10935 T 
V. parahaemolyticus LMG 2850 T 
•V. alginolyticus LMG 4409 T 
V. mytili LMG 19157 T 
Wpe/ag/usATCC25916 T 



m 



Wfort/sLMG21557 T 

V. cincinatiensis LMG 7891 T 
V. metschinilkovii LMG 11664 T 
-Kordatf/LMG13544 T 
— V. anguillarum LMG 4437 T 
V. diazotrophicus LMG 7893 T 

V. hispanicus LMG13240 T 
V. navarrensis LMG 15976 T 
Kvu/n/ffcusLMG13545 T 

V. ichthvoenteri LMG 19664 T 
— Vscophthalmi LMG 19158 T 
■ V aestuarianus LMG 7909 T 
^Wc/iagasi/LMG 21353 T 

E-V sptendidus LMG 19031 T 
V, tasmaniensis LMG 2001 2 T 
— V. lentus R3884 
V pomeroyi LMG 20537 T 

' K kanaloae LMG 20539' 

V. orientalis LMG 7897 T 
V. mediterranei LMG 11258 1 
K tapetis LMG 19706 T 

K pectenicida LMG 19642 T 



FIGURE 1 | Subtree incongruence test. Using the 58 vibrio taxon data set reported in Sawabe et al. (2007), subtree topology was compared between 8 gene 
(— 16S rRNA) (A) and 9 gene (+16S rRNA) phylogeny (B). The topology was also compared between 8 gene and only 16S rRNA gene phylogeny (C). 
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8-gene phylogeny corresponded to the clades that we previously 
proposed based on the 9-gene phylogeny (Figure IB). The results 
of the subtree incongruence test using the 58 vibrio taxa data 
showed that the inclusion of the 16S rRNA gene sequence is not a 
critical factor in optimizing the vibrio phylogeny on the basis of 
MLSA. 

THE LATEST VIBRIO PHYLOGENY BASED ON MULTILOCUS 
HOUSEKEEPING PROTEIN-CODING GENE SEQUENCES 

WGSs of key vibrios species that were resistant to the gene 
amplification, e.g., V. gazogenes, S. costicola, V. porteresiae, 
V. caribbenthicus, are now available. This result indicated that 
we could use the complete set of 8 housekeeping protein- 
coding gene sequences currently available from 86 described 
vibrio species and 10 genome-sequenced Vibrionaceae strains for 
the MLSA updating of the 8-gene phylogeny, on the basis of 
Splits Decomposition Analysis (SDA) (Bandelt and Dress, 1992; 
Huson and Bryant, 2005) (Figure 2) and a supertree recon- 
struction (Figure 3). On the basis of SDA, we could retain the 
14 distinct monophyletic clades that were previously defined, 
and we were able to further define 8 new clades: Damselae, 



Mediterranei, Pectenicida, Phosphoreum, Profundum, Porteresiae, 
Rosenbergii, and Rumoiensis (Figures 2, 3, Table 1). The robust- 
ness of these clades was high enough to propose their monophyly 
in the supertree reconstruction using three different molecu- 
lar phylogenetic analyses (Figure 3). Using an 8-housekeeping 
protein-coding gene analysis, most of clades shared >80.5% ANI 
and >92% AAI, and the highest ANI (98.3%) was observed in 
the sequence comparison between V. anguillarum and V. ordalii 
(Table 1). 

NEW CLADES 

Mediterranei consisted of three species: V. mediterranei, V. mar- 
itimus, and V. variabilis. The 8-gene ANI and AAI were 89.5- 
96.3% and 98.0-99.1%, respectively. The mol% G+C range of the 
clade members was 42-46.3 mol%. The genome-sequenced strain 
V. mediterranei AK1 showed 98.7% ANI and 99.9% AAI against 
the V. mediterranei type strain. Their known habitats are warm 
seawater and coral mucus. 

Porteresiae consisted of V. porteresiae and the newly described 
V. tritonius sp. nov. Detailed information for this new species 
is described in the section "The bacterial taxonomical remarks 



Pectenicida Scopthalmi Nereis Orientalis 



\ Coralliilyticus 



\Harveyi 



Gazogenes 




Mediterranei 

Nigripulchritudo 



/ Splendidus 



Salinivibrio-Grimontia-Enterovibrio Rosenbergii Phosphoreum 

Profundum Damselae 



Fischeri 



FIGURE 2 | Concatenated split network tree based on eight gene loci. 

The ftsZ, gapA, gyrB, mreB, pyrH, recA, and topA gene sequences from 96 
taxa were concatenated, and a tree was reconstructed using the SplitsTree4 



program. Clades indicated by a solid red line were the "new" clades 
proposed in this study. Clades indicated by a dotted red line or by a dotted 
black line are the clades "emended" and "un-changed," respectively. 
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-E. norvegicus LMG 19839' 
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W proteolytics LMG 3772' 
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V tenuis R3884 
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V. ichthyoenteri LMG 19664' 

V. nereis LMG 3895' 
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— K sinaloensls DSM 21326' 
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FIGURE 3 | Supertree reconstructed on the basis of the same data set drawn for Figure 2. 
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Table 1 | Newly proposed and emended clades by means of 8 gene MLSA for vibrios. 



Clade 


All described species included in newly 
proposed or newly involved species in 
the emended clade (in gothic) 


Total 
No. of 

species 


Genome 
sequenced 
strains involved 
in the related 
clade 


8 gene MLSA 
concatenated 
identity (%) 


8 gene-AAl 

(%) 


Habitat 






Damselae 


P damselae subsp. damselae and 
P. damselae subsp. piscicida 


Z 




yo.z 


QQ 7 


SsawatGr and fish 


Mediterranei 


V. mediterranei, V. maritimus, and 
V. variabilis 




AM 


oy.b— 9D.O 


QQ ft QQ 1 

sjo.u— yy. i 


Seawater and 
mucus of healthy 
coral 


Pectenicida 


V. caribbeanicus and V pectenicida 


2 




82.8 


96.1 


Marine sponge and 
diseased bivalve 
larvae 


Phosphoreum 


P. angustum, P iliopiscarius, P. leiognathi, 
and P. phosphoreum 


4 




87.5-95.8 


95.7-99.3 


Seawater, luminous 
organ, and fish 


Profundum 


R profundum, P. indicum, and P. lipolyticum 


3 




- 


- 


Deep sea 


Porteresiae 


V. porteresiae and V. tritonius sp. nov. 


2 




91.4 


97.2 


Wild rice and the 
gut of sea hare 


Rosengergii 


P. lutimaris and R rosenbergii 


2 








Seawater, coral, 
and tidal flat 


Rumoiensis 


V. litoralis and V. rumoiensis 


2 




91.5 


98.3 


Tidal flat sediment 
and drain of a fish 
processing factory 


ORPHAN CLADE 














Tapetis 


V. tapetis 


1 








Clam 


Proteolyticus 


V. proteolyticus 


1 








The intestine of 
isopod 


F6 




1 


F6 








F10 




1 


F10 










Anguillarum 


V. anguillarum, V. aestuarianus, and 
V. ordalii 


3 




87.0-98.3 


95.8-99.8 


Brackish, seawater, 
and fish 


Cholerae 


V. cincinnatiensis, V. cholerae, V. furnissii, 
V. fluvialis, V. metschnikovii, V. mimicus, 
and V. parilis 


7 




83.4-94.4 


93.5-99.8 


Brackish, seawater, 
and clinical specimen 


Coralliilyticus 


V. coralliilyticus and V. neptunis 


2 




94.2 


99.5 


Seawater, bivalve, 
and rotifer 


Diazotrophicus 


V. diazotrophicus and V. hispanicus 


2 




89.6 


97.0 


Brackish and 
seawater 


Fischeri 


Ali. fischeri, AIL logei, Ali. salmonicida, 
Ali. sifiae, and Ali. wodanis 


5 




88.7-93.5 


94.9-99.1 


Seawater, squid, and 
fish 


Gazogenes 


V. aerogenes, V. gazogenes, 
V. rhizosphaerae, and V. ruber 


4 




80.5-92.2 


92.0-99.0 


Estuary, salt marsh 
mud, and 

rhizosphere of wild 
rice 



(Continued) 
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Table 1 | Continued 



Clade 


All described species included in newly 
proposed or newly involved species in 
the emended clade (in gothic) 


Total 
No. of 
species 


Genome 
sequenced 
strains involved 
in the related 
clade 


8 gene MLS A 
concatenated 
identity (%) 


8 gene-AAl 

(%) 


Habitat 


Halioticoli 


V. breoganii, V. ezurae, V. gallicus, 

V. halioticoli, V. inusitatus, V. neonatus, 

and V. superstes 


7 




85.7-97.5 


95.3-99.6 


Gut of abalone and 
bivalve larvae 


Harveyi 


V. alginolyticus, V. azureus, V. campbellii, 
V. communis, V. harveyi, V. mytili, 
V. natriegens, V. parahaemolyticus, and 
V. rotiferianus 


9 


ATCC_BAA1116, 
EJY3, and Ex25 


86.7-96.0 


96.9-99.9 


Seawater, salt marsh 
mud, marine animal, 
and mucus of the 
coral 


Nereis 


V. nereis and V. xuii 


2 




89.9 


96.9 


Seawater and srimp 


Nigripulchritudo 


V. nigripulchritudo and V. penaeicida 


2 








Seawater and srimp 


Orientalis 


V. brasiliensis, V. hepatarius, V. orientalis, 
V. sinaloensis, and V. tubiashii 


5 




88.8-93.7 


97.1-98.9 


Brackish and 
seawater 


Scopthalmi 


V. ichthyoenteri, V. ponticus, and 

V. scophthalmi 


3 


N418 


87.8-94.4 


96.7-99.4 


Gut of flat fish, and 
marine animal 


Splendidus 


V. chagasii, V. crassostrea, 

V. cyclitrophicus, V. fortis, V. gigantis, 

V. kanaloaei, V. lentus, V. pelagius, 
V. pomeroyi, V. splendidus, and 
V. tasmaniensis 


12 


LGP32 


89.1-97.3 


96.9-99.9 


Seawater and marine 
animal 


Vulnificus 


V. navarrensis and V. vulnificus 


2 




86.9 


96.5 


Sewage, seawater, 
eel, and oyster 


Salinivibrio- 
Grimontia- 
Enterovibrio 
(Super clade) 


S. costicola subsp. costicola, G. hollisae, 
E. coralii, and E. norvegicus 


4 


AK16 and 1F-23 


0 81.7-89.4 


92.4-98.0 


Brine, human 
feceace, gut of 
turbot, and seawater 



*Bold indicated new and/or emended information. 



of V. tritonius sp. nov." below. These two species shared 91.4% 
ANI and 97.2% AAI. Two of the unique phenotypes in these 
species were an efficient H2 production and nitrogen fixation. 
While the genome sequences of these two species are highly con- 
served (unpublished data), they have distinct habitats (Table 1). 
The mol% G+C ranged from 44.2 to 45.5, and the DDH value of 
V. tritonius type strain against V. porteresiae type strain was 59% 
(Table 2). 

Pectenicida consisted of two species, V. caribbeanicus and 
V. pectenicida showing 82.8% ANI and 96.1% AAI. The reported 
habitats were tidal flats and diseased larvae, respectively (Table 1). 

Rumoiensis consisted of two species, V. litoralis and V. rumoien- 
sis showing 91.5% ANI and 98.3% AAI. These species were 
isolated from a tidal flat and sewage from a fishery product fac- 
tory, respectively (Table 1). The reported DDH value between 
V. litoralis and V. rumoiensis was below 7.4%. 

Damselae, Phosphoreum, Profundum, and Rosenbergii were the 
newly proposed clades that are included in Photobacterium spp. 
These four new clades are based on ANI (87.5-96.2% in range), 



AAI (95.7-99.3% in range), and branch separation according 
to the supertree analysis in comparison with those ranges and 
branch separations of other Vibrio clades. The Damselae clade 
consisted of two subspecies of P. damselae. The Rosenbergii clade 
consisted of P. lutimaris and P. rosenbergii. 

DEFINING ORPHAN CLADES 

These are the clades that are formed by only one species. V. tapetis 
and V. proteolyticus were not grouped with any other species 
(Figure 2 and Table 1). The recently proposed genomospecies F6 
and F10 also did not belong to any of the clades proposed in 
this analysis. Previously reported singletons, V. agarivorans, and 
V. pacinii were not included in this analysis due to the lack of some 
gene sequences. 

EMENDED CLADES 

We can find emendations in most of the clades previously 
defined (Figures 2, 3, and Table 1): (1) Cholerae (inclusion of 
V. parilis); (2) Fischeri (incl. AIL sifae); (3) Gazogenes (incl. 
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V. rhizosphaerae); (4) Halioticoli (incl. V. breoganii, and V. inusita- 
tus); (5) Harveyi (incl. V. azureus, and V. communis); (6) Orientalis 
(incl. V. sinaloensis); (7) Scophthalmi (incl. K ponticus); and (8) 
Splendidus (incl. K cyclitrophicus, V. gigantis, and V crassostrea) . 

We first included the complete set of the 8-housekeeping 
protein-coding gene sequences of Salinivibrio costicola subsp. cos- 
ticola in the MLSA for vibrio phylogeny, because its WGS data 
(ASAIO 1000001) are available. However, for the current analy- 
sis of the Salinivibrio/ Grimontia/Enterovibrio grouping, we could 
use only the data set including a single species of Salinivibrio, 2 
species of Enterovibrio, and a single species of Grimontia for the 
8-gene phylogeny. As both SDA and supertree analysis showed 
less robustness in these genera/species grouping, we decided to 
tentatively define them as the Salinivibrio-Grimontia-Enterovibrio 
(SGE) super-clade. 

DEFINING THE CLADE OF GENOME SEQUENCED STRAINS 

Vibrios are one of the most advanced groups in WGS analy- 
sis; currently more than 900 genomes are available in the public 
database (http://www.vibriobiology.net/). The MLSA of the 10 
genome sequenced strains revealed: (1) LGP32 and EX25 formed 
a robust cluster with V. tasmaniensis and V. alginolyticus, respec- 
tively; (2) N418 and EJY3 (Roh et al, 2012) were related to 
V. scopthalmi and V. natriegens, respectively; and (3) The orphan 
positions of genomospecies F6 and F10, and AK16 (Figures 2, 3, 
and Table 1). 

THE BACTERIAL TAX0N0MICAL REMARKS OF Vibrio tritonius sp. nov. 

On the basis of the 8-gene MLSA, the sea hare {Aplysia kuro- 
dai) isolates were highly likely to represent a new species within 
the family Vibrionaceae, more precisely within a new clade 
"Porteresiae." Four strains of V. tritonius formed a robust cluster 
within the Porteresiae clade on the basis of 4-gene sequence SDA 
(data not shown). To confirm the taxonomic status of the sea hare 
strains, a standard polyphasic taxonomy was conducted. 

The results of our phylogenic analyses based on the 16S 
rRNA gene sequence clearly showed that these strains belong 



to class Gammaproteobacteria, and more precisely to the fam- 
ily Vibrionaceae. The closest phylogenic neighbor of the four sea 
hare isolates was the V. furnissii-V. fluvialis cluster (Figure 4). 
V. porteresiae was not closely related, as shown by the 16S rRNA 
gene sequence phylogeny. Intra-species sequence similarities of 
the 16S rRNA gene among V. tritonius sp. nov. were above 99.5%. 
Four strains of V. tritonius sp. nov. showed 98.0-98.1% similarity, 
and 98.1% similarity toward V. furnissii (X76336) and V. fluvi- 
alis (X74703), respectively. Sequence similarities of V. tritonius sp. 
nov. to the other phylogenetic neighbors and to gas-producing 
vibrios were below 98%. The 16S rRNA gene sequence similarity 
between V. tritonius sp. nov. and "Allomonas enterica" AJ550855 
was 98.3%. The 16S rRNA gene sequences of V. fluvialis X74703 
and "A. enterica" AJ550855 were identical. 




V. tritonius sp. nov. JCM 16458 (GU95 1 700) 
V tritonius sp. nov. JCM 16459 (OU951701) 
| ■ V. tritonius sp. nov. JCM 16457 (GU951699) 
' — V. tritonius sp. nov. JCM 16456 T (GU951698) 
V. fluvialis ATCC 33809 T (X74703) 
V. furnissii ATCC 35016 T (X76336) 

V. brasiliensis LMG 20546 T (AJ316172) 
V. tubiashii ATCC 19109 T (X74725) 
V. hepatarius LMG 20362 T (AJ345063) 
V. mil LMG 2 1 346 T ( AJ3 1 6 1 8 1) 
V. nereis ATCC 259 1 7 T (X747 1 6) 

V. aemgenes ATCC 700797 T (AF124055) 

V. gazogenes ATCC 29988 T (X74705) 

100 I V. rhizosphaerae MSSRF3 T (D0847 1 23) 

I V. ruber JCM 11486 T (AF462458) 

V. porteresiae MSSRF30 1 (EF488079) 
V. cholerae CECT 5 14 T (X76337) 



FIGURE 4 | Unrooted phylogenetic tree on the basis of 16S rRNA gene 
sequences. Scale bar: 0.005 accumulated change per nucleotide. This 
figure combines the results of three analyses i.e., neighbor-joining, 
maximum parsimony, and maximum likelihood. The topology shown was 
obtained using neighbor-joining and 500 bootstrap replications. Percentages 
indicate the branches that were also obtained both in the maximum 
likelihood analysis (P < 0.01) and in the most parsimonious tree. 



Table 2 | DNA relatedness among Vibrio tritonius and the related vibiro species. 



Strain 




G+C content 
(moles %) 




% Reassociation with biotinylated DNA from: 




V. tritonius 
JCM 16456 T 


V. tritonius 
JCM 16457 


V. furnissii 
LMG 7910 T 


V. hepatarius 
LMG 20362 T 


V. tritonius 


JCM 16456 T 


44.2 


100 


92 


5 


5 


V. tritonius 


JCM 16457 


44.3 


74 


101 


NT 


23 


V. tritonius 


JCM 16458 


45.5 


76 


104 


12 


11 


V. tritonius 


JCM 16459 


45.2 


70 


100 


13 


14 


V. porteresiae 


MSSRF 30 T 


44.4 


59" 


NT 


NT 


NT 


V. furnissii 


LMG 7910 T 


52.0 


12 


15 


100 


14 


V. fluvialis 


LMG 7894 T 


49.3-50.6 


20 


NT 


NT 


NT 


V. hepatarius 


LMG 20362 T 


46.0 


10 


9 


3 


100 


V. nereis 


LMG 3895 T 


48.0 


6 


13 


4 


20 


V. tubiashii 


LMG 10936 T 


46.1 


8 


7 


7 


24 



"The reciprocal DDH value of V. tritonius JCM 16456 T against V. porteresiae MSSRF 30 T probe was 46%. 
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Mutual DDH experiments showed that the four strains of 
V. tritonius sp. nov., JCM 16456 T , JCM 16457, JCM 16458, and 
JCM 16459, were conspecific and clearly separated from their 
phylogenetic neighbors, e.g., V. porteresiae, V. fluvialis V.furnissii, 
V. tubiashii, V. hepatarius, and V. nereis (Table 2). The mol% G+C 
content was 44.8 ± 0.6, which was within the range of the genus 
Vibrio. 

DESCRIPTION OF Vibrio tritonius sp. nov 

Etymology of the newly describing Vibrio species was provided 
here: Vibrio tritonius (tri.to'ni.us. L. masc. adj. tritonius, named 
after Triton (a sea-god, son of Neptune and the nymph Salacia, 
referring to the habitat of the bacteria). 

Major phenotypic features of V. tritonius sp. nov. are shown 
in Table 3. The four sea hare strains have the major phenotypic 
features of the genus Vibrio (except for no growth on TCBS and 
gas production). These strains required salt for their growth, and 
they were motile, fermentative and oxidase positive. Apparent 
catalase activity was not observed. The four strains of V. trito- 
nius sp. nov. were phenotypically most similar to V. porteresiae, 
but they differed from V. porteresiae in four traits (catalase pro- 
duction, and the assimilation of D-mannose, y-aminobutyrate 
and pyruvate), out of 62 tested traits (Table 3). The four V. tri- 
tonius strains were sensitive to the vibrio-static agent 0/129 
(150 |ig). Positive assimilation of glucose, mannitol, gluconate, 
glucuronate, and xylose indicated the presence of three major car- 
bohydrate metabolic pathways, the Embden-Meyerhof, Entner- 
Doudoroff, and pentose-phosphate pathways, in V. tritonius sp. 
nov. Presence of the gene set for those three central metabolic 
pathways of carbohydrates was supported by our preliminary 
WGS analysis of V. tritonius JCM 16456 T (data not shown). 
Phenotypic traits differentiating V. tritonius sp. nov. from V. aero- 
genes, which shows a gas production phenotype, included nitrate 
reduction, amylase production, and arginine dihydrolase activity. 
Inability to grow on TCBS was a common trait of V. tritonius sp. 
nov. and V. porteresiae (Table 3). 

The other phenotypic traits were also described below. No 
swarming cells were observed. Gas production from glucose and 
mannitol occurred. Cells are curved rods, with rounded ends, 
are 0.7-0.9 |im in diameter and 2.6-2.7 |im in length when the 
organism is grown on ZoBell 2216E medium; the cells occur 
singly on the agar. No endospores or capsules are formed. 
Colonies on ZoBell 2216E agar medium are beige, circular, and 
smooth and convex with an entire edge. Sodium ions are essential 
for growth. The bacterium can grow in presence of 0.5% to 6% 
NaCl. The bacterium is a mesophilic chemoorganotroph which 
grows at temperatures between 15 and 40° C. Optimal growth is 
observed from 25 to 30°C. Growth occurs from pH 4.5 to pH 
9, and optimal growth is at pH 7.5-8.0. No growth occurs at 
45°C The bacterium is positive for acid production from glu- 
cose and mannitol; for nitrate reduction, acetoin production, and 
hydrolysis of gelatin, DNA and casein. The bacterium also can 
assimilate N-acetyl-D-glucosamine, cellobiose, D-fructose, mal- 
tose, D-mannitol, D-galactose, lactose, L-glutamate, L-proline, 
acetate, citrate, fumarate, DL-malate, pyruvate, and succinate. 
The bacterium is negative for catalase; indole production; argi- 
nine dihydrolase, lysine decarboxylase, ornithine decarboxylase, 



Table 3 | Phenotypic characteristics for distinguishing Vibrio tritonius 
from the related Vibrio species. 



Characteristics 


V. tritonius V. porteresiae V. furnissii 


V. fluvialis 




In =4) 


MSSRF 30 T 


LMG 7910 T 


LMG 7894 T 


Gas production 


+ 


+ 


+ 




from D-glucose 












37°C 


+ 


+ 


+ 


+ 


40°C 


+ 


+ 


+ 


+ 


Catalase 




+ 


+ 


+ 




Amylase 






+ 


+ 


Lipase 


- 


+w 


+ 


+ 


Acetoin production 


+ 


+ 






Arginine dihydrolase 






+ 


+ 


Indole production 






+ 


+ 


Growth on TCBS 






+(Y) 


+(Y) 


UTILIZATION OF 










D-Mannose 


d + (1) 




+ 


+ 


Melibiose 


d + (3) 








Lactose 


+ 


+ 






D-Glucuronate 


+ 


+ 




+ 


Trehalose 






+ 


+ 


Putrescine 






+ 




y— Aminobutyrate 


d + (3) 




+ 


+ 


Propionate 






+ 


+ 


Maltose 


+ 


+ 


+W 


+ 


D-Glucosamine 


+ 


+ 


+ 


+ 


D-Xylose 


+ 


+ 






Cellobiose 


+ 


+ 


+w 


+ 



Data were obtained under same culture condition in our laboratory. +, 
Positive; —, negative; w, weak reaction; d+, variable but type strain positive; 
d—, variable but type strain negative (numbers of positive strains are shown in 
parentheses). All species are motile, require Na for growth, fermentative of D- 
glucose, growth at 15-25° C, oxidase-positive, produce gelatinase, DNase, and 
caseinase, reduce nitrate, produce acid from D-glucose, assimilate N-acetyl- 
D-glucosamine, D-fructose, D-galactose, gluconate, D-glucose, D-marmitol, L- 
glutamate, L-proline, acetate, citrate, fumarate, DL-malate, pyruvate, succinate. 
All species are negative for pigmentation, swarming, luminescence, growth 
at 4 and 45° C, alginase production, agarase production, lysine decarboxylase, 
ornithine decarboxylase, requirement for organic growth factor, assimilation of 
D-sorbitol, aconitate, a-ketoglutarate, L-tyrosine, meso-erythritol. 

luminescence, and pigmentation; the requirement of organic 
growth factors; hydrolysis of agar, alginate, starch, and Tween 
80; and assimilation of D-glucosamine, D-sorbitol, aconitate, a- 
ketoglutarate, L-tyrosine, meso-erythritol, trehalose, putrescine, 
propionate, and D-glucosamine. The G+C content of DNA 
is 44.2-45.5 mol%. The type strain is JCM 16459 T = LMG 
2540 1 T = AM2 T . 

DISCUSSION 

Considerable biodiversity can be found within the family 
Vibrionaceae (Gomez-Gil et al., 2014), even after the first pro- 
posal of vibrio phylogeny and evolution was inferred on the basis 



Frontiers in Microbiology | Aquatic Microbiology 



December 2013 | Volume 4 | Article 414 | 10 



Sawabe et al. 



Updating Vibrio clades 



of MLSA in 2007 (Sawabe et al, 2007). More than 60 species of 
Vibrionaceae, with a surprising level of biodiversity, have been 
described since 2007. These include a marine invertebrate isolates 
such as coral associated vibrios (Chimetto et al, 2011; Gomez- 
Gil et al., 2014), introduction of nitrogen-fixing vibrios within 
an endophyte-like ecological niche (Rameshkumar et al, 2008), 
and an isolation of new vibrio species from the surface of cheese 
(Bleicher et al., 2010) have been reported. In addition to the 
increasing number of newly described vibrio species, many strains 
showing interesting ecophysiological features have been genome 
sequenced. However, taxonomic information appears to be insuf- 
ficient to push the elucidation of vibrio biodiversity and evolution 
forward. Such a rapid progress in the study of vibrio biodiver- 
sity, genomics and evolution prompted us to update the vibrio 
phylogeny on the basis of MLSA. We have retrieved the com- 
plete sets of 8 house-keeping protein-coding gene sequences for 
30 additional Vibrionaceae species including a newly described 
vibrio species, V. tritonius sp. nov., as well as for 10 as yet unnamed 
Vibrio/ Enterovibrio spp. The MLSA led us to propose eight 
new clades (Damselae, Mediterranei, Pectenicida, Phosphoreum, 
Profundum, Porteresiae, Rosenbergii, and Rumoiensis) in the fam- 
ily Vibrionaceae, in addition to those previously proposed in the 
report of Sawabe et al. (2007). In 2007, V. mediterranei, V. peteni- 
cida, V. rumoiensis, and P. rosenbergii were affiliated as singlet 
species, but they have now been grouped. Four orphan clades 
{Tapetis, Proteolytics, F6 and F10) were newly defined. More 
efforts are required to isolate the closest neighbors of these orphan 
species. Strains EJY3, EX25, N418, and LGP32 clustered robustly 
with V. natriegens, V. alginolyticus, V. scopthalmi, and V. tasman- 
iensis, but a further systematic survey is required to analyze the 
phylogenetic position of AK16 and 1F-230 (Figure 3). 

We are still facing a lack of Photobacterium spp. sequences 
to infer their precise evolutionary history. Among the 23 
described Photobacterium spp., we could include only half 
of them in this study. This situation has arisen mainly due 
to the "primer problems" in MLSA. Unfortunately, there are 
also limited numbers of WGS of Photobacterium spp. available 
in public databases. However, in this analysis, considering 



the results of SDA, supertree analysis, and the ANI and AAI 
similarity ranges in comparion to the other Vibrio spp. clades, 
we proposed four new clades for the Photobacterium spp.; (1) 
Damselae, (2) Phosphoreum, (3) Profundum, and (4) Rosenbergii. 
A Salinivibrio I Grimontial Enterovibrio super-clade is also 
proposed. 

In molecular phylogenetics, the use of minimum gene set is 
crucial to reduce time and cost, as well as to improve the accuracy, 
of analyses. This is of particular importance when identifying 
species and elucidating population structure and evolution in a 
super bacterial taxon such as the family Vibrionaceae, which has 
more than 140 species. The previous MLSA of 58 vibrio taxa 
(Sawabe et al., 2007) showed that 16S rRNA gene sequences have 
an extremely low species/strain discriminating power compared 
to the other genes tested. Therefore, before conducting the cur- 
rent vibrio MLSA, we evaluated whether the 16S rRNA gene 
data set could be eliminated from the MLSA. For this analysis, 
we developed a subtree incongruence test algorithm. The algo- 
rithm is a fast and reliable method for selecting subtrees that 
share the same topology or those that have different topologies 
but share the same species. The results of this analysis indicate 
that inclusion of 16S rRNA gene sequences is not necessary for 
reconstructing the vibrio phylogeny on the basis of MLSA. 

We have experienced the first case in which the molecu- 
lar phylogenies resulting from 16S rRNA gene sequences and 
from housekeeping gene sequences were largely incongruent 
in the species descriptions of V. porteresiae and V. tritonius 
(Figures 2, 4). For the affiliation of the clade of V. portere- 
siae, we used only four genes (the pyrH, recA, rpoA, and 16S 
rRNA genes), and we confirmed that V. porteresiae was affili- 
ated with the Cholerae clade (Rameshkumar et al, 2008).The 
16S rRNA gene sequence phylogeny revealed that V. tritonius 
sp. nov. was the most closely related to V. furnissi and V. ful- 
vialis with ca. 98% sequence similarity, and V. tritonius sp. 
nov. and V. porteresiae were distantly related in their phyloge- 
netic relationship (Figure 4). Less phylogenetic relatedness in 
the 16S rRNA gene sequence tree between both the species 
and the lack of housekeeping gene sequences of V. porteresiae 



Table 4 | FAME dominance (%) of Vibrio tritonius and the related species. 


Major fatty 






V. tritonius 




V. porteresiae 


V. fluvialis 


V. furnissii 


acids 


AM2 T 








MSSRF 30 T#1 


LMG 7894 T 


LMG 7910 T 




MA12 


MA17 


MA35 








Cl2:0 


3.5 


3.0 


3.0 


3.0 


6.8 


2.2 


2.1 


Cl4:0 


7.7 


6.5 


6.1 


5.7 


8.2 


3.5 


3.6 


Cl6:0 


24.6 


30.0 


32.8 


29.1 


18.8 


21.6 


19.2 


Cl8:0 


1.7 


0.8 


1.0 


0.9 


0.6 


2.4 


3.1 


Cl4 ; 1 (o7c 


2.0 


1.0 


0.8 


1.1 




tr 


0.5 


C 16:1 co7c 


28.5 


32.0 


33.7 


33.0 


27.5* 1 


36.8 


32.8 


C 18: ia>7c 


8.3 


17.8 


4.8 


19.1 


19.2 


26.9 


32.2 


11-methyl Ci8:iio7c 


7.8 


0.5 


2.8 


Tr 


0.82 


tr 


tr 


C 12:0 3-OH 


5.3 


4.6 


4.0 


4.6 


6.4 


3.2 


3.2 


Cl4:0 3-0 H 


2.0 


1.7 


1.4 


1.5 


6.8* 1 


1.1 


1.1 


Unknown ECL 20.454 


5.0 


ND 


5.6 


ND 




0.8 


ND 



* 1 Data from Rameshkumar et al. (2008). Values in Cj6 ju>7c and Cu:0 3-OH are from the summed percentage of feature 2 and 3, respectively. 
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prevented no direct comparison of V. porteresiae and V. trito- 
nius sp. nov. Fortunately, the whole genome nucleotide sequences 
of V. porteresieae and V. tritonius sp. nov. were determined in 
this study, and the first direct comparison of both species by 
MLSA and whole genome comparison was achieved. Surprisingly, 
the MLSA with 8-housekeeping genes phylogeny led to the 
conclusion that both Vibrio species, V. porteresiae and V. tri- 
tonius sp. nov., share a common ancestry and that they can 
be proposed as a new vibrio clade, "Porteresiae." Our pre- 
liminary genome comparison of both species also supported 
monophyly because a strong synteny was observed between both 
genomes (unpublished data). Incongruences between the 16S 
rRNA gene sequence tree and the MLSA tree were also observed 
in Mediterranei and Pectenicida clades (Lambert et al., 1998; 
Chimetto et al., 2011; Hoffmann et al., 2012) in this study. 
Therefore, we conclude that at present the "8-housekeeping-gene 
phylogeny" is the most powerful method for delineating vib- 
rio species description/biodiversity/population study/evolution, 
until alternative genome-based approaches are proposed. This 
analysis reduces the misidentification of ancestry clades of vibrios. 

The polyphasic taxonomic approach reveals that the Aplysia 
gut isolates, V. tritonius sp. nov., to be a novel Vibrio species show- 
ing gas producing ability (Figures 2-4, Tables 2, 3) and forming 
a robust clade, Porteresiae, with V. porteresiae. Production of 
gas during the fermentation of carbohydrates is not a prevalent 
property in the Vibrio genus (Shieh et al., 2000, 2003; Farmer 
et al., 2005; Kumar and Nair, 2007; Rameshkumar et al, 2008). 
V. aerogenes, V. furnissii, V. gazogenes, V. porteresiae, V. ruber, 
and V. rhizosphaerae are the species in which the gas production 
phenotype is retained in a stable way among 98 Vibrio species. 
For this reason, gas production is an atypical property in the 
genus Vibrio. The current phylogenetic network analysis using the 
8-gene MLSA confirmed that V. aerogenes, V. gazogenes, V. rhi- 
zosphaerae, and V. ruber form a robust clade, Gazogenes (Figure 2) 
(Kumar and Nair, 2007; Sawabe et al, 2007). V. furnissii belongs 
to the Cholerae clade. The newly described V. tritonius sp. nov. 
belongs to a new clade, Porteresiae. The gas compositions of 
Porteresiae clade species but also of the Gazogenes and Cholerae 
clades, were identical for H2 and CO2, but the H2 production 
efficiencies differ between these clades. The H2 production effi- 
ciency of V. tritonius sp. nov. and V. porteresiae was high and 
comparable to that of enterobacterial species, such as Escherichia 
coli, Salmonella, Enterobacter, and Klebsiella (Nakashimada et al., 
2002). Our preliminary genome comparisons suggests that V. tri- 
tonius sp. nov. and V. porteresiae contain very similar gene clusters 
responsible for H2 production and nitrogen fixation machinery 
(unpublished data). It would be intriguing to understand how 
those vibrios acquire and/or lose gas producing abilities from 
the standpoint of both evolutionary dynamics and metabolic 
diversity. 

In conclusion, 8-gene MLSA is a reliable tool for delineating 
a species and a monophyletic group or "clade." Using the cur- 
rent data set reported in this study, <98% of 8-gene-concatenated 
nucleotide sequence identity may allow us to define a species 
boundary. We also showed that WGSs overcome the limitations of 
gene-by-gene multilocus sequencing tasks found in the MLSA for 
Vibrionaceae (primer problem, Chromosome 2 gene inclusion). 



In fact, we could not include V. aerogenes, V. gazogenes, V. rhi- 
zosphaerae, and V. superstes for the previous 9-gene MLSA by the 
primer problems (Sawabe et al., 2007), but the successful WGS 
for these species and the inclusion of 8 housekeeping protein gene 
sequences retrieved from the genome sequences can provide the 
better picture of current vibrio molecular phylogeny. More efforts 
to sequence the individual housekeeping genes and WGS of all 
remaining species in the family Vibrionaceae not involved in this 
study allow the possibility of elucidating the ultimate clade struc- 
ture of these Vibrios. The ultimate phylogenetic trees allow us 
to provide the ideal phylogenetic backbone to elucidate the evo- 
lutional history, genome dynamics, and plasticity in the family 
Vibrionaceae. 
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